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Abstract 



' In this letter we will use higher-order supersymmetric quantum mechanics to obtain several fam- 

ilies of complex solutions g{x;a,b) of the Painleve IV equation with real parameters a,b. We shall 
also study the algebraic structure, the eigenfunctions and the energy spectra of the corresponding 
non-hermitian Hamiltonians. 

Keywords: quantum mechanics, non-linear differential equations, Painleve equations, complex 
' potentials with real spectra 

» ' 

1 Introduction 

' Since its birth, supersymmetric quantum meclianics (SUSY QM) catalyzed the study of exactly solvable 

^ , Hamiltonians and gave a new insight into the algebraic structure characterizing these systems. Histori- 

0^ ' cally, the essence of SUSY QM was developed first as Darboux transformation in mathematical physics [T] 

and as factorization method in quantum mechanics [2^ . 

On the other hand, there has been an increasing interest in the study of non-linear phenomena, 
which in many cases leads to the analysis of Painleve equations Oil]. Although these were discovered 
from strictly mathematical considerations, nowadays they are widely used to describe several physical 
phenomena [5]. In particular, the Painleve IV equation {Piv) is relevant in fluid mechanics, non-linear 
optics and quantum gravity |^. 

As it has been shown, there is a natural connection between quantum systems described by second- 
order polynomial Heisenberg algebras (PHA), whose Hamiltonians have the standard Schrodinger form 
^ \ and their differential ladder operators are of third order, and solutions g{x\ a, b) of Piv HUZIIH]- Moreover, 

^ ■ these algebras can be realized by the fc-th order SUSY partners Hf, of the harmonic oscillator Hamil- 

- - -' tonian Hq, which leads to a simple method for generating solutions of P/y, since the SUSY technique 

provides explicit expressions for the extremal states of Hk and connecting formulae relating them with 
the corresponding Pjy solutions [9HT2]. It is worth to note, however, that the need to avoid singularities 
in the new potential Vk{x) and the requirement for the Hamiltonian Hk to be hermitian lead to some 
restrictions jl3) : (i) first of all, the relevant transformation function has to be real, which implies that 
the associated factorization energy is real; (ii) as a consequence, the spectrum of Hk consists of two inde- 
pendent physical ladders, an infinite one departing from Eq = 1/2 (the ground state energy of Hq) plus 
a finite one with k equidistant levels, all of which have to be placed below Eq. Regarding Piv, these two 
restrictions imply that non-singular real solutions g{x; a, b) can be obtained for certain real parameters 
a, b. 

From the point of view of spectral design, it would be important to overcome restriction (ii) so that 
some (or all) steps of the finite ladder could be placed above Eq- In this way we would be able to 
manipulate not just the lowest part of the spectrum (as we did previously [Tni[ITl[T3]), but also the 
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excited state levels, which would endow us with new tools for spectral design. In this article we are going 
to show that this can be achieved if one relaxes restriction (i) as well, which will force us to use complex 
transformation functions (see ^M) and will lead to the generation of complex solutions to P/y. 

This letter is organized as follows: in Section 2 we shall present the general framework of SUSY QM 
and PHA. In the next Section we will generate the complex solutions to P/v, specifically, we shall analyze 
the domain of the parameter space (a, b) for which we obtain real or complex solutions; then, in Section 
4 we will study the eigenfunctions and the energy spectra of the non-hermitian Hamiltonians. We shall 
present our conclusions in Section 5. 



2 General framework of SUSY QM and PHA 



In the k-tYi order SUSY QM one typically starts from a given solvable Hamiltonian 
and generates a chain of standard (first-order) intertwining relations [15H17| 



^ V2 



l,...,/c. 



Hence, the following equations must be satisfied 



(1) 

(2) 
(3) 
(4) 



(5) 

(6) 



We are interested in the final Riccati solution afc(a;,efe), which turns out to be determined, either by fc 
solutions Q!i(a;, Cj) of the initial Riccati equation 



^■^[x,tj) + a{{x,t.j) = 2[yo(a;) - e^], j = l,...,k, 



or by k solutions Uj of the associated Schrodinger equation 



Hou-j 



u'- + Vq{x)uj = ejUj, j = 1 



, . . . , fc. 



(7) 



(8) 



with ai{x,ej) — u'j/uj. 

Thus, there is a pair of k-th order operators interwining the initial Hq and the final Hamiltonians iJfe, 
namely, 

(9) 



where 



HkB. - B. Ho, HqB, - B, Hk, 



B+ = A+...At, B^^A^...A-,. 



(10) 



The normalized eigenfunctions of Hk, associated to the eigenvalues Em are proportional to the 
action of B^ onto the corresponding ones of Hq (■0„, n = 0, 1, . . . ). Moreover, there are k additional 

eigenstates V'cj'' associated to the eigenvalues Ej (j = I, . . . , k), which are simultaneously annihilated by 
B^ . Their corresponding explicit expressions are given by |13[|18) : 



(fc) _ 



[(P„-ei )...(£;„-€,.)] 1/2' 

W{ui, . . .,Uj^i,Uj+i, ...,Uk) 
W{ui,...,Uk) 



E„ 



(11) 
(12) 
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Let us note that, in this forniahsm the obvious restriction < i?o = 1/2 naturally arises since we want 
to avoid singularities in Vfc(x). 

On the other hand, a m-th order PHA is a deformation of the Heisenberg-Weyl algebra of kind [TT1IT8] : 



[i/,L±] =±i±, 

[L-, L+] = Q„, + i(i/ + 1) - Qrn+l{H) = P„,{H), 

771 

Q,n+i{H) = L+L- =l[{H-£,), 



(13) 
(14) 

(15) 



where Qm+i{x) is a (m + l)-th order polynomial in x, which implies that Pmix) is a polynomial of order 
m in X and Si are the zeros of Qm+i{H), which correspond to the energies associated to the extremal 
states of H. 

Now, let us take a look at the differential representation of the second-order PHA (m = 2). Suppose 
that is a third-order differential ladder operator, chosen by simplicity as: 

(16) 
(17) 



(18) 



L+ 


= lX 


L 


+ 

2 1 




Li 


1 








V2 


dx 


1 




1 






-f h(x) 


~ 2 





These operators satisfy the following intertwining relationships: 

HL\ = LX{Hs_ + 1), HaL2 = L2H, 
^ \H,L+]=L+, 



(19) 
(20) 



where 77a is an auxiliary Schrodinger Hamiltonian. Using the standard equations for the first and second- 
order SUSY QM and decoupling the resulting system gives rise to 



f ^x + g 
h 



where 



5 = ^ + 0.9' + W + 2[x' ^a)g+-. 
Note that this is the Painleve IV equation {Piv) with parameters 

a = E2+S:i-2Ei-l, b = -2{£2 - £3)^ 



2 , 5' 

+2 


9' 


2 


g' 


9^ 


Y ~ 2" 


+ y 


3 3 
+ 2^ + 





xg 



(21) 
(22) 

(23) 
(24) 

(25) 



Hence, if the three quantities £i are real, we will obtain real parameters a, b for the corresponding P/y. 

Let us note that the potential of equation (j23l) contains a harmonic oscillator term. In addition, three 
terms can be identified {—g'/2 + g^ /2 + xg) which in general lead to an anharmonicity in the potential 
and are completely determined by the solution to P/y. As a consequence, one could say that the solution 
g to P/y is the main responsible for the spectral differences which the Hamiltonian H could have with 
respect to the harmonic oscillator (compare with [TU]). 
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3 Complex solutions to Pjv with real parameters 

It is well known that the first-order SUSY partner Hamiltonians of the harmonic oscillator are naturally 
described by second-order PHA, which are connected with P/y, as we have shown in the previous section. 
Furthermore, there is a theorem stating the conditions for the hermitian higher-order SUSY partners 
Hamiltonians of the harmonic oscillator to have this kind of algebras (see [l3j). The main requirement is 
that the k Schrodinger seed solutions have to be connected in the way 

Uj = [a-y-^ui, (26) 
e, 61 - (j - 1), j = l,...,fc, (27) 

where a~ is the standard annihilation operator of Hq so that the only free seed ui has to be a real solution 
of Eq. dHI without zeros, associated to a real factorization energy ei such that ei < — 1/2. 

In this work we intend to overcome this restriction, although if we use the formalism as in |13j with 
ei > Eq, we would obtain only singular SUSY transformations. In order to avoid this we will instead 
employ complex SUSY transformations. The simplest way to implement them is to use a complex linear 
combination of the two standard linearly independent real solutions which, up to an unessential factor, 
leads to the following complex solutions depending on a complex constant X + in (A, k G M) [14]: 



u{x; e) 



72 



^ , l-2e 1 2\ /^ ■ X ^ /^3-2e 3 2 
1^1 ( -^'2' J + ^^(^ + 1-^1 [^^'2' 



(28) 



where iFi is the confluent hypergeometric (Kummer) function. The known results for the real case [20] 
are obtained by making k — and expressing A as 

^ = ^''Ti^^y ^^^^ 

with e M. 

Hence, through this formalism we will obtain the fc-th order SUSY partner potential Vk{x) of the 

harmonic oscillator and the corresponding Pjy solution g{x; ei), both of which are complex, in the way 

Vk{x) = Y-{HWiu,,...,Uk)]y', (30) 

g{x;e,) = -x-{Hijs,{x)]y. (31) 

Note that the extremal states of Hk and their corresponding energies are given by 

W{ui,. . . ,Uk-i) , . , , 

1 °^ ~fF7 T' '^^i = efc = £1 - A:- 1 , (32) 

Ty(Mi, . . . ,Ufe) 

i,^,^Bte'^"l\ E2^\, (33) 

^^£3 oc B+a+ui, ^3 = ei + 1. (34) 



Recall that all the Uj satisfy Eq. (|26|) and ui corresponds to the general solution given in Eq. ((28)) . 

For fc = 1 , the first-order SUSY transformation and Eq. (I3ip lead to what is known as one-parameter 
solutions to P/y, due to the restrictions imposed by Eq. (|25p onto the two parameters a, h of Piv which 
makes them both depend on ei |21j . For this reason, this family of solutions cannot be found in any 
point of the parameter space (a, 6), but only in the subspace defined by the curve {(a(ei), 6(ei)) , ei G M} 
consistent with Eqs. ([^5)) . Then, by increasing the order of the SUSY transformation to an arbitrary 
integer fc, we will expand this subspace to obtain k different families of one-parameter solutions. This 
procedure is analogous to iterated auto-Backlund transformations [22]. Also note that by making cyclic 
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Figure 1: Parameter space (a, &) of the Piv solutions. The curves represent the solution subspace for 
non-singular real or complex (solid curves) and only complex (dashed curves) solutions. 



permutations of the indices of the three energies 8i and the corresponding extremal states of Eqs. p2j|34p . 
we expand the solution families to three different sets, defined by 



2 

2' 

02 = 2ei - fc, 62 = -2fc^ (36) 



ai = -ei+2fc-^, 6i = -2(ei + i) , (35) 



3 , „ / . 1^ ^ 



a^^-ei-k--, 63 = -2 (^ei + , (37) 

where we have added an index corresponding to the extremal state (Eqs. (I32j|34p ) taken to build up the 
Piv solution in Eq. (j3ip . The first pair, ai, &i, can provide non-singular real or complex solutions, while 
the second and third ones can give just non-singular complex solutions, due to singularities in the real 
case. A part of the non-singular solution subspace for both real and complex cases is shown in Fig.[T] One 
can check that those points which belong to two different sets have associated the same Piv solutions. 

In turn, let us analyze some of the Pjy solutions obtained by this method. The real solutions arise by 
taking k = 0, and expressing A as in Eq. with ei < 1/2. They can be classified into three relevant 
solution hierarchies, namely, confluent hypergeometric, complementary error and rational hierarchies. Let 
us note that the same set of real solutions to Piv can be obtained through inverse scattering techniques [5] 
(compare the solutions of j2l| with those of [l3]). In Fig. [21 three real solutions to Piv are presented, 
which belong to the complementary error hierarchy. 

Next, we study the complex solutions subspace, i.e. we allow now that ei > 1/2. The real and 
imaginary parts of the complex solutions g{x] a, b) for two particular choices of real parameters a, 6, 
which belong to different solution sets, are plotted in Fig. |3l 

Note that, in general, ip£^ ^ V x S M , i.e., the solutions g{x; a, b) are not singular. Moreover, both 
real and imaginary parts have an asymptotic null behaviour (g — >■ as \x\ — >■ 00). This property becomes 
evident in Fig. [31 as well as in the parametric plot of the real and imaginary parts of g(x] a, b) of Fig. [31 



4 Non-hermitian Hamiltonians 

Let us analizc the Hamiltonian Hk obtained by the complex SUSY transformation. Note that the real 
case, which leads to hermitian Hamiltonians, has been studied previously fTSpS] . allowing to obtain some 
criteria related to the structure of the associated energy spectrum Sp(iJfc), the number of zeros of the 
eigenfunctions of H^, and the way in which they are connected by the third-order ladder operators L^. 
This action is in agreement with the fact that Sp{Hk) consists of an infinite ladder plus a finite one: there 
are two extremal states (both annihilated by L^) from which the two ladders start, one associated to Ck 
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Figure 2: Some real solutions to Piv, corresponding to oi = 1, 6i = (fc = 1, ei = ~l/2, ly — 0.7) 
(solid curve), Oi = 4, 5i = —2 (fc = 2, ei = —3/2, i/ = 0.5) (dashed curve), and ai = 7, 6i = —8 (A; = 3, 
ei = —1/2, v = 0.3) (dotted curve). 




Figure 3: Real (solid curve) and imaginary (dashed curve) parts of some complex solutions to P/y. The 
upper plot corresponds to a2 = 12, 62 = — 8 (fc = 2, ei = 7, A = K = 1) and the lower one to as = —5, 

&3 = -8 (A; = 1, ei = 5/2, A = k = 1). 




-4-20 2 4 

Re[g(x;a,b)] 



Figure 4: Parametric plot of the real and imaginary parts of g{x;a,h) for ai = —6, 61 — —2 (fc = 1, 
ei = 5/2, A = 1, K = 5) and |a;| < 10. For bigger values of x, the curve approaches the origin in both 
sides. 
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Figure 5: Spectrum of the SUSY partner Hamihonians (right) and Hk (left) for ei > 1/2, ei 7^ Ej] 
Sp(-fffc) still contains one finite and one infinite ladders. The dark bars represent the original and mapped 
eigenstates of iJo and Hk , while the light ones the k new levels of Hk introduced by the fc-th order SUSY 
transformation. All of them have associated square-integrable eigenfunctions. 



and the other one to Eq = 1/2; since the ladder starting from ends at ei, the eigenfunction associated 
to ei is annihilated by . The actions of onto any other eigenstate of Hk are non-null, and connect 
only the eigenstates belonging to the same ladder. 

As far as we know, complex SUSY transformations with real factorization energies were used for the 
first time by Andrianov et al. to obtain non-hermitian Hamiltonians with real spectra [14 . These topics 
have been of great interest in the context of both parity-time (PT) symmetric Hamiltonians (see Bender 
et al. [24]) and pseudo-hermitian Hamiltonians (see Mostazafadeh et al. [25]). Next, we will examine the 
structure of non-hermitian SUSY generated Hamiltonians Hk- 

First of all, the new Hamiltonians necessarily have complex eigenfunctions, although the associated 
eigenvalues are still real. In previous works, the factorization energy associated to the real transformation 
function ui was bounded, ei < Eq — 1/2. In this paper we are using complex transformation functions 
to be able to overcome this restriction and yet obtain non-singular solutions. This naturally leads to 
complex solutions to Pjv generated through factorization energies which could be placed now above Eq . 
The resulting spectra for the non-hermitian Hamiltonians Hk obey the same criteria as the real case, 
namely, they are composed of an infinite ladder plus a finite one, which now could be placed, either 
fully or partially, above Eq- The eigenfunctions associated to the energy levels of the original harmonic 
oscillator are given by Eq. pT|) and the ones associated to the new energy levels by Eq. (fT^ . all of them 
square-integrable. A diagram of the described spectrum is shown in Fig. [5l 

The extremal states of the SUSY generated Hamiltonian Hk are given by Eqs. (|32ll34p . These are 
non-singular complex eigenfunctions of Hk and, from their asymptotic behaviour, we conclude that those 
given by Eqs. p2|33p are square-integrable. Note that in this case the oscillatory theorem does not hold 
anymore, neither for the real nor for the imaginary parts, although a related node structure emerges. 
The absolute value and the real and imaginary parts of tpg-^ (x) for two particular cases are shown in Fig. 

m 

On the other hand, complex transformations for ei = Ej are worth of a detailed study, namely, when 
the factorization energy ei belongs to the spectrum of the original harmonic oscillator Hamiltonian. For 
instance, let us consider a first-order SUSY transformation with ei = Ej and ui given by Eq. (|28p. i.e., 
ui is a complex linear combination of the eigenfunction ijjj of Hq and the other linearly independent 
solution of the Schrodinger equation. It is straightforward to see that the action of the ladder operator 
L~ = Afa~A^ is given by 

El, L-iA+^ji)^A+^i^i, (38) 
E,, L-{AttlJj)^Q, (39) 
Eo, L-{A+ ^0)^0, (40) 
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Figure 6: Plot of the absolute value, the real, and the imaginary parts (solid, dashed and dotted lines, 
respectively) of the eigenfunction ips^ {x) given by Eq. ([5^ for the values k = 2, ei = —1, A = 1, k = 1/2 
(up) and fc = 2, ei = 4, A = 1, K = 6 (down). 



where I ^ j, I ^ 0, the shown energies correspond to the departure state, and we have used that 
A^ipj (X l/ui. For = A^a^A^ we have 

El, L+(A+^i)«V^i+i, (41) 
E„ L+{A+yj,) = 0, (42) 

which does not match with the established criteria for the non-singular real and complex cases with 
ei 7^ Ej since now it turns out that: 

E,+i, L-(A+^,+i)cxA+^, cx — ^0, (43) 

Ej^i, L+{A+i;j_i) (X A+^Pj (X — ^0. (44) 

Ul 

The resulting Hamiltonian is isospectral to the harmonic oscillator but with a special algebraic structure 
because now one state (the one associated to Ej) is connected just in one way with the adjacent ones 
(associated to Ej ±1). A diagram representing this structure is shown in Fig. [T] We are currently 
studying the fc-th order case and expect to find the new criteria which will be valid for these special 
transformations. 



5 Conclusions 

In this letter, based on PHA and higher-order SUSY QM, we have introduced a method to obtain real and 
complex solutions g{x;a,b) of the P/y with real parameters a, 6. We have studied the properties of the 
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Ho Hi 




Figure 7: Spectra of the harmonic oscillator Hamiltonian Hq and of its first-order SUSY partner Hi when 
we use the factorization energy ei = Ej e Sp(iJo)- The level Ej of Hi is connected with its adjacent ones 
just in one way. 

resulting solutions, including the analysis of the subspace of the parameter space (a, h) were non-singular 
real or complex solutions can be found. In addition, we have analyzed the algebras, the eigenfunctions 
and the spectra of the non-hermitian SUSY generated Hamiltonians. 

Further investigation on the description of the analytic structure of the complex solutions and on 
the way in which the ladder operators map between different eigenstates and their corresponding 
energy levels is needed. Besides, we are looking for extensions of this technique to obtain Pjy solutions 
associated to complex parameters a, h. 
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